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Abstract. We discuss unusual aspects of symmetry that can happen due to entropic 
effects in the context of multi-scalar field theories at finite temperature. We present 
their consequences, in special, for the case of nonrelativistic models of hard core 
spheres. We show that for nonrelativistic models phenomena like inverse symmetry 
breaking and symmetry non-restoration cannot take place, but a reentrant phase at 
high temperatures is shown to be possible for some region of parameters. We then 
develop a model of interest in studies of Bose-Einstein condensation in dilute atomic 
gases and discuss about its phase transition patterns. In this application to a Bose- 
Einstein condensation model, however, no reentrant phases are found. 
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1. Introduction 

One of the most interesting aspects concerning the studies of multi-field models at finite 
temperature is the possibility of emergence of a much richer phase diagram than one 
would usually find in one field type of models. The possibility that unusual symmetry 
patterns could emerge in those models, for some specific region of parameters, has 
attracted considerable attention in the literature (see, for instance, [TJ |2j and references 
therein and [Sj for a short review). 

For most of the standard physical systems we know in nature, we have a good sense 
of how symmetries seem to change as the temperature is changed. Typically, the larger 
is the temperature the larger is the symmetry exhibited by the system and vice-versa. 
Examples of this behaviour are expected to happen in the particle physics models, like in 
the electroweak phase transition and possibly also in Grand-Unified models. The same is 
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expected in the much lower energy systems, like those of condensed matter. For example, 
the phase transition in ferromagnets, superconductors, Bose-Einstein condensation of 
atomic gases, etc, just to name a few systems. In all these examples we always go from 
an ordered (less symmetric) phase below some critical temperature of phase transition to 
a disordered (more symmetrical) phase above the critical temperature, or the opposite, 
if the temperature is decreased from a high temperature (or symmetry restored) state 

s 

However, the above symmetry aspects seem not to be the rule. In fact we are 
also becoming increasingly aware that entropic effects in multi-field models may show 
other patterns of symmetry breaking and restoration that are less usual. For example, 
many condensed matter systems, like spin glasses, compounds known as the manganites, 
liquid crystals and many others, commonly show phenomena like reentrant phases of 
lower symmetries at higher temperatures and, therefore, they can exhibit unusual phase 
diagrams that we would otherwise not expect. Many of these systems have recently 
been reviewed in As concerned to quantum field theory models, the possibility of 
other phase transition patterns was also shown to be possible in the context of multi- 
scalar field theories at finite temperatures pQ. These models show the possibility of a 
symmetry that is not broken at low temperatures, getting broken at high temperatures 
(what is called an inverse symmetry breaking). Other case that seems possible is a 
symmetry that is broken at lower temperatures, not getting restored at all as we go to 
higher temperatures (what is called a symmetry non-restoration). The problem of how 
a symmetry broken or restored phase may emerge as a reentrant phase in the system, 
like it is seen in many low energy condensed matter systems, was recently analyzed in 
j2j in the context of a coupled nonrelativistic model of two scalar fields. 

The plan of this paper is as follows. In the next section we briefly review the results 
obtained in [2] for the case of a nonrelativistic model of two scalar fields with overall 
symmetry U(l) x £7(1) and show that it admits reentrant phases for some region of 
parameters. One possible physical realization of this kind of model is for example in 
the description of a coupled two species dilute atomic gas system. This is a kind of 
system in atomic physics that has been of great interest recently concerning studies 
(both theoretical and experimental) of Bose-Einstein condensation. In Sec. 3 we offer 
a quantum field theory description for this problem. We then analyze the possibility 
of emergence of reentrant phases in the quantum field formulation for these systems. 
This would be a novel symmetry behaviour that could be of great interest, given its 
possible implementation in the laboratory. Our results allow us to conclude, at the level 
of our approximations, on the non appearance of such reentrant phases in these kind of 
coupled dilute atomic gases systems. 

2. Reentrant Phases in Nonrelativistic Multi-Scalar Field Models 

We start our discussion by considering the following nonrelativistic Lagrangian density 
model, of two (complex) scalar fields $ and with global symmetry U(l) x U(l), 
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1 T-l2\ ,Tr ,, VTV*VTY ^* ATr*iTA2 



V 2m^ / 3! 

(1) 

The model (0) can be thought as coming from the nonrelativistic limit of a corresponding 
relativistic counterpart, as shown explicitly in In the interaction parameters g$ 
and g^ describe two-body self-interaction terms, as commonly considered for dilute and 
cold (low energy) systems of particles 0, in which case only binary type interactions, 
i.e., hard core type of interactions of the form shown in are relevant, g is the 
cross-coupling between the two fields, which we consider here as a quadratic type of 
interaction. The mass parameters m$ and are the masses for the fields (or particles). 
The one-body type of interactions, of magnitude and fi^, can either represent the 
effect of external potentials (for example a magnetic field) on the system, internal energy 
terms (like the internal molecular energy relative to free atoms in which case the fields 
in the Lagrangian would be related to molecular dimers), an explicit gap of energy in 
the system (like in superconductors), or just chemical potentials added to the action in 
the grand-canonical formulation to enforce finite density (or fixed number of particles) 
for both $ and The latter will be the case for our application of Q to the coupled 
atomic gas problem in Sec. 3. Here we will just consider the \i\ parameters as constant 
one-body parameters added to our model such that the possible symmetry breaking 
patterns, depending on the sign of fii (i = $, can easily be determined from the 
potential term in (£[]). This is just what is done in spontaneous symmetry breaking 
studies performed on the relativistic analogous models. Therefore, for /ij > 0, we 
have an initially symmetry restored phase in both $ and \1/ directions, while fii < 
corresponds (at zero temperature) to symmetry broken phases for both $ and 

We require the model to be overall bounded from below, which then gives 
the constraint condition on the two-body interaction terms, g^ > 0, g§ > and 
g^g<s> > 9g 2 . This is the same condition imposed on the analogous relativistic problem 
PUI2|. This boundness constraint will be observed in all our results below. Non-trivial 
phase transitions can emerge for negative values of the cross-coupling g, which is allowed 
by the above boundness constraint. This was shown in diverse instances to be the case 
in the relativistic analogous models [TJ |3] . 

In the analysis below, we will also restrict, for simplicity, to an initial symmetry 
restored phase (at zero temperature) for both fields, \ii > and leave the symmetry 
broken case for the Bose-Einstein condensation problem studied in Sec. 3. The phase 
structure of the model is then determined by the sign of the temperature dependent 
one-body terms, /ij(T) = + Sj, where Sj is the field temperature dependent self- 
energy. We look for an intermediate (reentrant) phase at some interval of temperature 
and parameters in which the symmetry in one of the field directions is broken. This 
analysis can in principle be carried out within perturbation theory, as described in 



Symmetry Aspects in Nonrelativistic Multi-Scalar Field Models and Application 



4 



a companion paper (3 a , which indeed shows the possibility of appearance of reentrant 
phases for a system of hard-core particles described by ([TJ. However, the results in [21 E] 
also shows that as the (perturbative) temperature corrections are considered for the two- 
body terms, g<s>(T), gy(T), these effective couplings run to negative values above some 
temperature T neg . This then violates the initial condition of boundness for the potential 
when considering the model in equilibrium in a thermal bath with temperature T > T neg . 
At the leading order perturbative calculation and assuming /!$ = fi^ = fi and m<$> = 

= m, T neg ~ min (l27r v //i/(2m 3 )^/(5 5 | + 9g 2 ), 12vryVV(2m 3 )W(%l + 9# 2 )). 
This is reminiscent of the breakdown of perturbation theory in quantum field theory at 
finite temperature, which is well known in relativistic models (see e.g. |Zj and references 
in there). Nonperturbative methods are then called for a proper interpretation of the 
results and to confirm that the appearance of reentrant phases in our model is not just 
an artifact of perturbation theory. 

The problem of the self-couplings running to negative values at high temperatures 
can be solved, e.g., by resumming all leading order bubble corrections to the couplings. 
This is naturally done in the context of the renormalization group, by solving the flow 
equations for all couplings and parameters of the model. A simpler and equivalent 
approach was also shown in [2j, where this resummation is also accomplished by 
solving a set of self-consistent homogeneous linear equations for all effective couplings, 
<?$(T), gm(T) and g(T), and the result of these equations feeded back in the equations 
for the effective one-body terms, fii(T). We refer the interested reader to |2J for the 
details and we give here only the main results of this approach. 

We consider, for illustrative purposes, the parameters (at T = 0): g<$> = 2 x 
10~ 15 eV~ 2 ,g = — 10~ 16 eV~ 2 , m$ ~ = 1 GeV and yU$ = /i^ = IneV. These values 
of couplings and masses could for instance be representative of some dilute Bose gas 
atom or molecule (see next section). The temperature and the tree-level value of g^ are 
then changed and we look for regions of symmetry broken phase in the \l/ field direction 
(for the values of parameters considered, it easy to show that the symmetry remains 
always restored in the $ direction |3J. For all values of parameters and temperature 
considered we check the boundness condition extended for the effective couplings 
(temperature dependent), as obtained by the flow equations described previously. That 
is, g*(T)g*(T) > 9g 2 (T), or that the ratio R NR (T) = g^(T)g^{T)/[9g 2 (T)} > 1. The 
resulting phase diagram, as a function of g$ at T = and the (log of the) temperature, 
is shown in Fig. 1. 

Fig. 1 shows clearly the possibility of reentrant phases in the system, through an 
inverse symmetry breaking, in \1/ direction. For instance, for gq,(T = 0) = 10~ 16 eV~ 2 , 
we find a reentrant symmetry broken phase starting at the temperature T ( f 1 ^ B ' ) ~ 
3.4 x 10~ 4 eV (or ~ 4 K) and ending at T C ( ^ R) ~ 1.4 x i(T 2 eV (or ~ 161 K), through 
symmetry restoration. In this region of temperature and parameters, i? NR (T) > 1 
and the (effective) potential is still bounded from below. In the $ direction there 
is no symmetry breaking or reentrant phases at any temperature for the parameters 
considered. 
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Figure 1. The phase diagram of the system in terms of and temperature, for the 
parameters considered in the text. The dark gray region denotes a reentrant phase 
with symmetry breaking (SB) in the 'J direction, /i$(T) < 0. The region below it, in 
light gray, is the unstable region, i? NR (T) < 1 and above it is a symmetry restored 
(SR) phase, /Lt*(T) > 1. Temperature is given in units of eV. 



3. Application to a Coupled Two-Species Dilute and Homogeneous Atomic 
Bose Gas Model 

Let us now consider the case of model (P) as describing two coupled Bose gases of 
fixed densities p$ and p^, respectively. The model could then be describing a system 
composed by a mixture of coupled atomic gases, like the ones recently produced [S], 
with same chemical element in two different hyperfine states, or even two different 
mono-atomic Bose gases in the homogeneous case Here /i<j> and p^ are then explicitly 
chemical potentials added in the grand-canonical formalism to ensure the fixed densities 
for each Bose atom gas. Here we start describing the system in the broken phase in both 
<3> and \1/ directions. Therefore, /i$ — > — p$ and /i* — > — /i^ in ((T} and these chemical 
potentials are taken as positive quantities and with their values determined by the usual 
thermodynamic relation, in terms of the pressure P(T, p<j,, /i^,), 

&P(T, <9P(T, /%) 
= ^ , P* ~ , {*) 

where the pressure is defined as the negative of the effective potential computed at its 
minima (which is the thermodynamic free energy of the system), 

P = P(T,//$,^) = -V eS (T, 0o, Vo) , , , , , (3) 

4>0=<Pm,1p0=1pm 

where m and ipm are the values of 0o and ipo that extremizes (corresponding to a 
minimum of) the effective potential, 
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dV eS (T,4> ,i/;o) 



d(j) 







4>0 = 4>m,4>0=lpn 



0. (4) 



The effective potential follows from (|TJ) by expanding the fields around the vacuum 
expectation values, (<&) = (po/y/2 and {$) = ijj Q /y/2, and it is evaluated in the one-loop 
approximation in a standard computation of quantum field theory at finite temperature 
(for the one- field case, see for instance JO]). At the tree- level, <p m and ip m are given in 
terms of the minima of the potential in (JTJ), 

,5"*//$ — Sg/i-is 2 „5"$A** — 3(7//$ 



in 



6- 



■0,1 



6- 



(5) 



3*0* - 9g 2 m g^g^ - 9g 2 

At finite temperature, the equations for m and ip m are given by analogous expressions to 
(0), but in terms of the effective chemical potentials instead, /i$ and p,^, that are defined 
by the solution of the self-consistent equations, /!$ = //$ — E^ and = /i^ — E^, 
given in terms of the $ and \l/ field self-energies E. 

The explicit expression for the pressure at finite temperature, obtained from 
the effective potential as described above and that follows from some length but 
straightforward calculation, is given by [TT] 



where 



and 



P(T, /t$, /ity) 



2 (0*0* - V) 



g® i/4 - e! 



d 3 q 

/3J (2tt) 3 



In (l - e _/3ji+ ) + In (l 



— e 



-/3A_ 



.4 



± 



■4>o = i>n 

<t>0 = 4>-n 



4>0 = V>m ' 
00 = 0m 



Hq/Gq, — H$G$) + 4:g 2 ip^ n (f) 2 n G<s,G^ 



1 1/2 



(6) 



(7) 



1 ,2 



2m$ 4 



2ma 



(8) 



.f/j and Gj denote the Higgs and Goldstone modes, respectively, for each field in the 
broken phase. 

From (J2J) and after some algebra to eliminate the dependence of these expressions 
on the chemical potentials, we find the expressions relating the total densities p# and 
p$ with (the condensate densities) ip m and <f> m , as given by [TT] 



it 1 

P *"-2 



(2^)3 



(l + 2n^ 



djjLq, 



(H-2nj 



(9) 
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and 



d 3 q 

(2tt) 3 



9Al / x dA„ i 

s^( 1 + 2 '^) + ^( 1 + 2 "' 



(10) 



where n^ ± = l/[exp(/L4±) — 1] and the partial derivatives of A± with respect to /i$ and 
fiq, are defined by 



dA A 
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+ G^ 




AA_ 


dA- 


(h^, + G^j 




aA_ 



H^Gi, — HgjGa, ) — 2g ^) m ^) m G^ 



H^Gy — H<s>G<i>) + 4g 2 4>l l ip^ n Ga ? G^ 



Hq,Gq, — H^Go}) + 2g 2 (J) 2 m il)'^ n G$ 



^9 2 $J) 2 m Gi>G^, 



H^G^j — H$G$ 
H^iG^j — HqGa} 



2g^ m r m G* 



Hq,Gq, - Hq,Gq, ) + ^(^J^Gi.G^ 



HtyGqi — H$G$) + 2g 2 (f) 2 n tl) 2 n G$ 



H^Gy — Hs>G<s>) + Ag 2 ^ 2 ^ 2 ^^^ 



(11) 



The coupled equations (jHJ) and (fTU|) give completely the phase diagram for the 
condensates ipm and m as a function of the temperature and the densities. A qualitative 
analysis of the phase structure is also possible to be deduced already at this level 
from equations ©, (fTU|) and (fTT]) . Note that for g = 0, from we obtain that 
A 2 + (g = 0) = H^G®, A 2 _(g = 0) = H^G^, and we obtain the Bogoliubov spectrum 
for each field in the uncoupled case. Also, (jHI) and (fTUJ) decouples and we obtain as a 
result, for example for p$, 



P "-^ + 2 



dh 



_31_ + 



</«< 



(27T) 



q 2 / q 2 i 

2m$ \ 2m$ 



(/*< 



1 + 



,(12) 



- 1J 

with analogous equation for p^. Taking (|T^|) at the critical point, T = T Cj $, we have 
that 4> m {T = T c $) = 0, since the condensate density at T c vanishes and (|T2^) gives p$ = 
/c/ 3 q/(27r) 3 [exp(q 2 /(2m$T Ci $)) - l]" 1 , or inverting it, T Ci<& = 27r/m$ [p$/C(3/2)] 2/3 , 
where C(3/2) ~ 2.612. This is the standard result for the critical temperature of an 
homogeneous ideal Bose gas. This result emerges because of the level of approximation 
we are considering. It is only modified by corrections due to the self-interactions through 
nonperturbative methods and it requires at least second order corrections in the self- 
energy (see, for instance, ^2] and references therein). 

Note also that, at the level of approximation we are considering, from the equations 
0, ©, (fTU|) and (JUJ), if any of the fields go above the transition point (either <p m = 
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or ip m = 0) the two equations (JHJ) and (fTUj) also decouple, becoming independent of each 
other, since the cross-coupling term in (J7J) and (fTTj) always appears multiplying both <ft m 
and ip m . As a result, no reentrant phase at high temperatures seems to be possible here. 
A computation performed in the restored phase case (similar to the one done in Sec. 
2) also seems to confirm this result. This comes about as a consequence of the strong 
temperature dependence introduced by the chemical potentials through the relation @ 
in both the broken and symmetric cases. In the broken (BEC) phase, it is also seen that 
the system exhibits a small dependence on the cross-coupling term (at one-loop order 
it is even insensitive to the sign of g). A throughout analysis of the phase structure 
coming from the coupled set of equations (JOJ) and (fTUJ) . including higher order terms, 
will be presented elsewhere 
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